Kinetic theory of QED plasmas in a strong 
electromagnetic field 
II. The mean-field description 



A. HolP' 1 V.G. Morozov b < 2 G. R6pke a ' 3 

a Physics Department, University of Rostock, Universitdtsplatz 3, D-18051 Rostock, 

Germany 

b Moscow State Institute of Radioengineering, Electronics, and Automation, 
117454 Vernadsky Prospect 78, Moscow, Russia 



Abstract 

Starting from a general relativistic kinetic equation, a self-consistent mean-field 
equation for fermions is derived within a covariant density matrix approach of QED 
plasmas in strong external fields. A Schrodinger picture formulation on space-like 
hyperplanes is applied. The evolution of the distribution function is described by 
the one-particle gauge-invariant 4x4 Wigner matrix, which is decomposed in spinor 
space. A coupled system of equations for the corresponding Wigner components is 
obtained. The polarization current is expressed in terms of the Wigner function. 
Charge conservation is obeyed. In the quasi-classical limit for the Wigner compo- 
nents a relativistic Vlasov equation is obtained, which is presented in an invariant, 
i.e. hyperplane independent, form. 

Key words: relativistic kinetic theory; QED plasma; hyperplane formalism; 
mean-field approximation 



1 Introduction 



In part one of this article [1] we have developed a covariant density matrix 
approach to kinetic theory of QED plasmas, making use of the relativistic 
hyperplane formalism [2-4] in the Schrodinger picture. In what follows the 
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paper [1] will be referred to as I and equations from this paper will be labeled 
by (I. . .), where ". . ." denotes the equation number. In the present paper we 
aim to derive quantum mean-field kinetic equations for the fermionic subsys- 
tem starting from a general relativistic covariant equation discussed in paper 
I. 

Section 2 briefly recalls notations and definitions from paper I needed in fur- 
ther considerations. Section 3 is devoted to the derivation of a mean-field 
kinetic equation for the gauge-invariant fermionic Wigner function in the hy- 
perplane formalism. We next use the spinor decomposition of the Wigner 
matrix and obtain a set of coupled covariant equations describing kinetic pro- 
cesses in different channels. The charge conservation is shown to be fulfilled. 
In Section 4 we discuss the quasi-classical limit in the hyperplane formalism. 
We derive mean-field kinetic equations for the distribution functions of par- 
ticles and antiparticles, and show that these equations can be represented 
in a fully covariant form. Section 5 will conclude and give a short outlook. 
In Appendix A we show the relation between our approach and the existing 
mean-field theories of QED [5,6]. Finally, Appendix B gives expressions for 
the matrix commutation and anticommutation relations which are necessary 
for the spinor decomposition of the kinetic equation. 

Except for the quasi-classical limit, we use the system of units with h = c = 1. 
The signature of the metric tensor is (+, — , — , — ). 



2 Basic definitions 

We consider a quantum plasma of charged fermions interacting through the 
electromagnetic (EM) field. For simplicity, we will take these fermions to be 
electrons and positrons, but the inclusion other fermions (say, protons) as 
additional Dirac fields is not a particular problem. The system is assumed to 
be subjected to a prescribed external EM field which is not necessarily weak. 

In paper I we defined the one-particle density matrix for fermions as the 

average 

Paa> (X ± ,X' ± ] T) = (p aa , ( X± , x' ± )) T = Tl {p aa , (x ± , X ± ) g(n, t)} , (2.1) 

where g(n, r) is the nonequilibrium statistical operator of the system on a 
hyperplane cr characterized by a unit time-like normal four-vector n M and a 
scalar parameter r = x ■ n which may be interpreted as an "invariant time" . 
The fermionic density operator p is given in the Schrodinger picture on the 
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hyperplane by 



= -l$ a M,h«)]> (2-2) 



a, a' being spinor indices of the Dirac field operators. The transverse compo- 
nents (which are space-like) of x — and other four-vectors V = {V^} are 
defined with respect to the normal n through the decomposition 

x P = n i* T + x P, V" = n"Vf, + Vf, (2.3) 

where 

V. = n ■ V, V? = A» U V U , A"„ = - n\. (2.4) 



Our further analysis rests heavily on the basic "equal-time" anticommuta- 
tion relations for the Dirac field operators on hyperplanes [cf. Eqs. (1.3.24) 
and (1.3.25)] 



(2.5) 



where 



»(xj = /^e^^p.,) (2.6) 



is the three-dimensional delta function on a hyperplane a and the matrix 
7n(n) is defined through the following decomposition of the Dirac matrices: 

(2.7) 



As discussed in paper I, the self-consistent mean-field approximation for the 
fermionic subsystem can be introduced only when the EM field variables are 
separated into the macroscopic condensate mode and the photon degrees of 
freedom. This we have shown in paper I by means of a time-dependent unitary 
transformation of the statistical operator and the operators of the EM field 
[see Eqs. (1.4.1) and (1.4.3)]. After this procedure, the effective Hamiltonian 
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describing the fermionic subsystem in the mean-field approximation can be 
taken in the form 

Hl{n) = J da -4 {-'-^M V.+m^-.+J daj^xJA^xJ, (2.8) 



where the symbol : O : shows the normal ordering in operators, and the space- 
like derivatives V M =V M — V M are defined by the relations 

= ^ + v, = a;^ = a;A (2 . 9) 



_L 



We will use the notation G(r, x ± ) = G{jit + x ± ) for any function G(x) on the 
hyperplane a n T furthermore. 

The first term in Eq. (2.8) is the Hamiltonian of free fermions, while the second 
term describes their interaction with the total mean EM field in the system, 
A 1 . We have shown in paper I that the total field tensor = d^A u — d v A 1 
satisfies Maxwell equations 

d^(x)=f(x)+f ext (x). (2.10) 
Here j" xt (x) is a prescribed external current and 

nx) = (j"(x ± )y (2.ii) 

is the mean polarization current. For the electron-positron plasma, the current 
operator is (e < 0) 

]"(x ± ) = e^(x ± h^(x ± ):. (2.12) 



If protons are treated as a dynamical subsystem, the corresponding term must 
be included into the current operator. 

As outlined in paper I, the total Hamiltonian of the system contains, in addi- 
tion to (2.8), the Hamiltonian of free photons, H EM (n), and the term Hl nt (n), 
describing the interaction between fermions and photons. In the mean-field ap- 
proximation the interaction term is neglected. The derivation of the fermionic 
kinetic equation requires to calculate commutators of the dynamical fermion 
operators ip and ip with the Hamiltonian (see below). The free photon con- 
tribution H EM (n), consisting completely of the photon dynamical operators 
d u A^, will not contribute to these commutators and can therefore be omitted. 
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This implies that the dynamics of the EM field in the mean-field approxima- 
tion is completely governed by the Maxwell equations (2.10). Nevertheless, at 
the end of the paper we shall discuss some non-trivial relations between the 
mean-field description of the fermionic subsystem and the photon kinetics in 
QED plasmas. 

To complete the list of definitions, we write down the expression for the gauge- 
invariant "one-time" Wigner function [9] on the hyperplane a n T . The Wigner 
function is expressed in terms of the one-particle density matrix (2.1) by 

W aa ,(x ± ,p ± ;r) = Jd%e^5(yn) 

x exp{ieA(x ± + ±y ± ,x ± - ly ± ;r)}p aa , (x ± + ±y ± ,x ± - \y ± ;rj , (2.13) 

where 

x ± 

A(x ± ,x' ± ;t) = J A^(r,R ± )dR» 

x' ± 

1 

= Jds « - xf) A^{t, x' ± + s(x ± - x'J) (2.14) 



is the gauge function. Our immediate task will be to derive a mean-field kinetic 
equation for W. 



3 Mean-field kinetic equations 

3. 1 Kinetic equation for the one-particle density matrix 

We start with the mean-field kinetic equation for the density matrix (2.1). 
Taking (2.8) as the effective Hamiltonian for the fermionic subsystem, we 
have 

The commutator in the right-hand side is easily calculated by using the iden- 
tity 
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[Paa'( x ±, x '±)i-My±)i j b'(y'±)-} = (i\Xb s3 ( x ± - y±) Pv a >(y±, x '±) 

- (7n) 6 v s3 ( x '± - y'±) Pab( x ±, y±), ( 3 - 2 ) 

which follows from the anticommutation relations (2.5). After some algebra 
we find that Eq. (3.1) can be written in matrix notation p = [p aa ,\ as 

-f v p(x ± ,x' ± ;r)_ 
+ ("^V M + eA ±fl (r, xj) S^p(x ± , x' ± ; r) 
+ (tV^ + eA^x'j) P (x ± ,x' ± ;t)S» 

-ie(A^(r,x ± ) - A l{ (r, x'j) p(x ± , x' ± ; r), (3.3) 

where 

S" = {P"n v , a^= l -[Y,Y}- (3-4) 

Eq. (3.3) defines the mean-field dynamics of the fermionic one-particle density 
matrix for time-like translations with respect to the space-like plane <J nr . 
In the special Lorentz frame where n M = (1, 0, 0, 0) (which sometimes will 
be referred to as the "instant frame"), we have r = t and therefore Eq. (3.3) 
describes the time evolution of the one-particle density matrix. It is interesting 
to note that, within the mean-field description, there is no need to know the 
explicit form of the nonequilibrium statistical operator g(n, r). The mean-field 
kinetic equation follows directly from the equation of motion for the density 
operator p with the effective Hamiltonian (2.8). 

3.2 Kinetic equation for the Wigner function 

Applying the Wigner transformation (2.13) to Eq. (3.3), we obtain in matrix 
notation W = [W aa ,] 

D T W = -im[ lv W] - i D XM [S», W] - P, {^, W} , (3.5) 
where we have introduced the operators 

1/2 

D T = ^--e J dsnT IJa ,(T,x ± -isV p )VZ, (3.6) 

-1/2 



—p(x ± ,x' ± ;T) = -im 
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1/2 

D ±M = V M -e J ds^ u (r,x ± - iS V p )v;, (3.7) 



-1/2 
1/2 

P»=P^-*e J sds^ u (r,x ± - iS V p )v;, (3.8) 

-1/2 

and the transverse gradient in the momentum space: = A' 1 " d/dp u ± . The 
transverse part of the total field tensor is defined as 

jrv» = - V v A*i. (3.9) 



The virtue of Eq. (3.5) is its compact and covariant form. It should be em- 
phasized, however, that the Wigner function, the matrices 7^ S^, and the 
operators (3.6) - (3.8) are defined with respect to the family of hyperplanes 
o nr characterized by the normal n M . The fact that n M is an arbitrary time-like 
unit vector reflects Lorentz covariance of Eq. (3.5). To see this, we note that 
the normal vectors of the same plane in different frames, n M and n' M , are re- 
lated by a Lorentz transformation (boost) = h} x v n v which is equivalent to 
the transformation of space-time coordinates x' M = A^ u x u . In the new Lorentz 
frame the invariant time parameter r has the same value, since r = n-x = n'-x'. 
For any given n, there exists the special "instant frame" where n M = (1, 0, 0, 0) 
and, consequently, r = t. In Appendix A we show that in this frame Eq. (3.5) 
reduces to the kinetic equation derived by Bialynicki-Birula et al. [5] on the 
basis of a different approach. 

Despite its apparently simple structure, Eq. (3.5) is a rather complicated ma- 
trix equation. To obtain a deeper physical insight into processes described by 
this equation, it is convenient to expand the Wigner function in a complete 
basis in spinor space 

W = \ (IW + 7/t W + 75 W (P) + 7 5 7 M W ( % + <v W") . (3.10) 

Here I is the unit matrix, and W M , W^, VV^, are the scalar, vector, 
pseudo-scalar, axial-vector and tensor coefficient function of the Wigner ma- 
trix W respectively. By the trace rules in spinor space it can easily be verified 
that the coefficient functions can be expressed as 



W = tr(W), (3.11) 
W M = tr ( 7 M W) , (3.12) 
W (P) =tr( 75 W), (3.13) 
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Wfa = tr ( 7 ^ 75 W) , (3.14) 
W = -tr (cr^W) , (3.15) 

where the symbol "tr" means the trace over spinor indices. Now a straight- 
forward algebra allows to derive a coupled set of equations for the coefficient 
functions from the kinetic equation (3.5) by using Eqs. (3.11) - (3.15) 



D T W = 2(n a P /J -n /J P a )W* (3.16) 

D T W = -(n^D ±a - n a D?)w a - 2e% x n a P^ A) - AmW^n a) (3.17) 

D T W (P) = 2imn a W? A) + 2ie af3Xg n a P^yV Xs : (3.18) 

D T W ( % = -2e% A ^ p/3 ^ A + 2tmn^W (P) - (n»D ±a - n a D») W? A) , (3.19) 
D r W^ = (raT" - nT") W - m{n^W u - n v W) + ie" ua(3 n a P p W {P) 

+ (n»D ±a - n a D?)W a - (n»D ±a - n a D»)w» a , (3.20) 

where e^ al3 is the Levi-Civita tensor. In Appendix B the necessary algebra 
including commutator and anticommutator relations is shortly surveyed. 

The tensor structure of Eqs. (3.16) - (3.20) becomes more clear if we split 
these equations into longitudinal and transverse components with respect to 
the hyperplane cr . The vector and axial- vector functions are decomposed 
according to 

W = n"W„ + Wf , Wf A) = n*W m + V^ (A) , (3.21) 

where 



(3.22) 



The tensor function can be written as 

= (WW - U u n») + Wf, 
W = n a W» a , VVf" = A^A^W^. (3.23) 
Then we arrive at the following set of equations: 

D T W = -4P a W°, (3.24) 
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D r W„ = 




(3.25) 


D T W? = 




(3.26) 


D r W (P) = 


= 2iroW ||(A) +2ie^n°P'Vl^, 


(3.27) 


D rW||(A) 

II v j v 


= 2tmW (P) - D ±a W? {A) , 


(3.28) 


D T WT(A) 


= -2e% x n a P^ + D»W }l{A) , 


(3.29) 


D T W = 


-P^W + mVVf - D ±Q >Vf tQ , 


(3.30) 


D T wr = 




(3.31) 



This representation for the mean-field kinetic equation will prove to be par- 
ticularly convenient in the quasi-classical limit. 



3.3 The charge conservation 



In order to describe the picture consistently the polarization current (2.11), 
which enters the Maxwell equations (2.10), must be expressed in terms of the 
Wigner function. As shown in paper I [see (1.5.9)], the polarization current 
can be written in the form 

r(x)=ej ^5(p-n)W^x ± ,p ± ;r = x-n), (3.32) 



where W M is the vector component (3.12) of the Wigner function. Let us 
prove that the fundamental charge conservation law d^j^ = is satisfied in 
our theory. 

It is convenient to define for any function G(x ± ,p ± ; r) the transformation 

G(x) = J -^-^5(p-n)G(x ± ,p ± ;r = x-n). (3.33) 



Then Eq. (3.32) takes a compact form 

f( x ) = eW(x). (3.34) 

It follows easily from Eq. (3.33) that 



d»G(x)=nJ^)+V fl G(x). (3.35) 
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Note also that for functions G which go to zero as \p ± \ — > oo we have 

D~6 = n* i d l P, D>[G = WG, WG = pjG. (3.36) 

Applying the transformation (3.33) to Eq. (3.17), we obtain 

- 2e e% x n a (p?Wfo) - 4emW^n,. (3.37) 
On the other hand, we may write 

V = ny^-" + v^-". 

Combining this with Eq. (3.37), we see that 9^ = 0. 

One can follow a similar procedure to derive other conservation laws and bal- 
ance equations for local quantities like the mass current, the spin density, the 
magnetic moment density, and the angular momentum density. The advantage 
of the hyperplane formalism over the previous approaches to the mean-field 
QED kinetic theory [5,6] is that all the balance equations and conservation 
laws will have a manifestly covariant form. 



4 The quasi-classical limit 

4-1 The local-field approximation 

In view of practical applications of the theory, it is of interest to study kinetic 
processes in QED plasmas depending on a slowly varying external EM field. 
To consider this case, we insert the constants H and c into Eq. (3.5) and the 
operators (3.6) - (3.8). Then we obtain the kinetic equation 

ITfiC 1 ' 

V T W = -— [ lv W\ - - D ±fl [5", W\ - - P„ {5", W} (4.1) 

and the following expressions for the operators with the corresponding gradient 
expansions: 

C -l/2 
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t - \ ^ ^ + ^ (V • V p ) 2 n"^ + . . . , (4.2) 

1/2 

V M - ^ y JF ±/i;/ (r, x x - is/iVj,) V p 

C -i/ 2 

V, " ^ ^ V£ + ^ ( V • V p ) 2 n"^ V* + . . . , (4.3) 
ieh f . _ / ■ - — ^ —71/ 



1/2 

f| (V • V p ) + ... (4.4) 



^ 12c 

The condition that the terms containing the field derivatives be small reads 

A B < ^emi (4-5) 

where A B is the average de Broglie wave length for fermions and l EM is the 
characteristic length for variations of the EM field in the system. For laser 
induced plasmas the latter quantity is roughly equal to the wave length of the 
external laser field. We will refer to the condition (4.5) as the local approxi- 
mation, from which the operators (4.2) - (4.4) can be concluded to be 



D r = J^- V^V;, (4.6) 



D ±M = V M - -JF^ , (4.7) 
P„=P^- (4-8) 



In the local-field approximation, Eqs. (3.24) - (3.31) (with inserted h and c) 
become 



D T W = -^ P±I M», (4.9) 
D T W„ = - D^VVT, (4.10) 

D T Wf = DfW,, - r^nVx^) - -=-«", (4.11) 



2mr 2 

D T W (P) = i— +^^a^>^, (4.12) 
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D r W ||(A) = i— W {P) - D ± ^ (A) , (4.13) 

D T WT (A) = -| e % A n>^W ± A + D£W II(A) , (4.14) 

| rji f 

V T W = --p^W + —W»- D ±a Wr, (4.15) 

ft ft 

D T Wf = ^ e^ Q/3 n aP±/3 W (P) - D£W + U^W. (4.16) 



It should be emphasized that the local approximation in the operators D r , Dx, 
and P M does not necessarily implies that all quantum effects are neglected. It 
can easily be seen from Eqs. (4.9) - (4.16) that some components of the Wigner 
function show non-analytic behavior in the limit H — > describing quantum 
phenomena like pair production in strong fields [7,8]. This aspect of the "one- 
time" mean-field theory of QED was discussed by Bialynicki-Birula et al. [5] 
in their study of the Dirac vacuum in strong external fields. For QED plasmas 
the situation is somewhat similar to the problem of the Dirac vacuum, but, 
generally speaking, a self-consistent description of quantum effect in plasmas 
involves the photon kinetics. This point will be detailed in a special article. 

4-2 Quasi- classical kinetic equations for fermions 

We now want to obtain the quasi-classical limit of Eqs. (4.9) - (4.16) for 
those components of the Wigner function which determine the polarization 
current (3.32). 

We first notice that Eq. (4.11) allows to eliminate U.^ in the other equations. 
In particular, we have 

which is to be inserted into Eq. (4.9). Then we can observe that, in the quasi- 
classical limit (h—>0), Eq. (4.15) leads to the relation 

if 

W^ = — W. (4.18) 

mc 

Thus Eqs. (4.9) and (4.10) give a closed set of quasi-classical equations for W 
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and W„: 



D r W - ^ P±M D r (K W) + -L K D ±M W„ = 0, 
D r W„ + -L D ±M (KW) = 0. 



With (4.6) and (4.7), it is easy to verify that in the quasi-classical limit 

Px,D/ = - V^X, D ± X = 0- ( 4 -20) 

Using these relations, a little algebra shows that Eqs. (4.19) can be written in 
a more symmetric form 

„ / , > x (4-21) 



c ^ \ mc £ 



where we have introduced the dispersion relation for fermions on the hyper- 
plane 



e(pj = c^/m 2 c 2 - p\ (4.22) 
and the transverse four-velocity 



c 2 



< = — rf. (4.23) 



We now define the distribution functions for electrons (w) and positrons (w) 
on the hyperplane a n T : 



w(x ± ,p ± ; t) = U ^ W(x ± , P± ; r) + W„ (x ± ,p ± ; r) \ , 



w(x x ,p ± ; t) = \ W(x ± , -p ± ; r) - W (a; x , -p x ; r) }> . 

2 



(4.24) 



These functions satisfy independent kinetic equations which follow from Eqs. (4.21). 
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In the expanded form, we have 



Yr + - V A w - - (n^ u + - jr\ V> = 0, 



These equations are in fact nothing more than a generalization of relativistic 
Vlasov equations to the case that the distribution functions for particles and 
antiparticles are defined on arbitrary hyperplanes a n T . In the special "instant 
frame", where = (1,0,0,0) and = (ct,r), Eqs. (4.25) take the well- 
known form [see Appendix A] 



dw _ 

— — h v ■ vw + e 

at 



dw 
~di 



+ v ■ Vw — e 



£ + -{v x B) 



£ + -(« x B) 

c 



— = 0, 
dp 



dw 
dp 



(4.26) 



0, 



where v = c 2 p/t(p) is the velocity vector, and e(p) = cyp 2 + m 2 c 2 is the 
relativistic energy in the "instant frame". 

The kinetic equations (4.25) describe the evolution of the distribution func- 
tions with respect to the time-like variable r. Note, however, that w and w 
can also be regarded as functions of the space-time point x and the four- vector 
p ± , according to 

w(x, p ± ) = w(x ± ,p ± ; t = x ■ n), w(x,p ± ) = w(x ± ,p ± ; r = x ■ n). (4.27) 



This interpretation of the distribution functions allows to put Eqs. (4.25) into 
a more compact form. First we note that 

which follows directly from the definition of J r ±fiU , Eq. (3.9). We next consider 
the relations 



d 



u^d tl w(x,p A ) =u>*V ll w(x,p ± ) + (u-n) \—w(x ± ,p ± ;T 



dr 
d 



u^d /x w(x : p ± ) = u^Vpivfapj) + (u-n) y—w(x ± ,p 



(4.28) 
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where the time-like unit vector is defined as 



= U + . (4.29) 

mc 2 \ c I 



Elimination of the r-derivatives between Eqs. (4.25) and (4.28) and some 
rearrangement leads to the equations 

(4.30) 

«" (^ + ^J>(x)V^)«;(x,p x ) = 0. 

One can verify, e.g., by going to the "instant frame", that (4.29) is the four- 
velocity of a particle with the four-momentum p^ satisfying the mass-shell 
condition p 2 = m 2 c 2 . Equations (4.30) may thus be interpreted as the evolution 
equations for the distribution functions with respect to the invariant proper 
time. 

In order to guarantee a self-consistent description of the plasma in the quasi- 
classical approximation, the polarization current (3.32) must be expressed in 
terms of the distribution functions w and w. To do this, we recall the quasi- 
classical result (4.18) and write 

if 

W = n"W,, + — W. (4.31) 

11 mc 

Elimination of W N and W with the aid of Eqs. (4.24) gives 

2 

TTIC 

WfopJ = ^—^{u^p ± )w(x,p ± )-u^-p ± )w(x,-p ± )}, (4.32) 

so that the polarization current (3.32) takes the form (with the inserted 
Planck's constant) 

/ct^p TflC 2 r ] 

(2-Khf S( " P ' ^ 7{p) ^ P± ) [ w ( x iP±) -w{x,p ± )\. (4.33) 



By using Eqs. (4.30), it can easily be verified that the above expression for 
the current is consistent with the conservation law d„j^ = 0. 
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4-3 The invariant quasi- classical distribution function for fermions 



It is interesting that the polarization current (4.33) can be rewritten in a 
form where the four-vector n does not appear. First we notice that the delta 
function 5(p ■ n) in Eq. (4.33) may be replaced by 5 (p ■ n — e(p ± )/c) because 
other functions in the integrand do not depend on p n — p • n. Then, according 
to the identity 



/ 



d A p 



5(p-n-e(p ± )/c) 
e(p±) 



(■■■) = - c J d 4 P 5(p 2 -m 2 c 2 ) (•••), (4.34) 

p°>0 



we can rewrite Eq. (4.33) as 

d 4 p 



j^(x) = 2emc 



p°>0 



5(p 2 — m 2 c 2 )u' J '(p ± ) w(x,p ± ) — w(x,p ± ) . (4.35) 



Finally, from Eqs. (4.23) and (4.29) follows 



u^(p ± ) = [p ■ n - e(p ± )/c]. 

mc mc 



(4.36) 



With the mass-shell constraint p 2 = m 2 c 2 we have 
Eq. (4.35) becomes 



f(x) = 2e 



/ 



d 4 p 
{2TTh) : 



f(x,p) - f(x,p) 



p^/mc, so that 



(4.37) 



after introducing the mass-shell distribution functions for particles and an- 
tiparticles: 



f(x,p) = Q(p°) 5{p 2 — m 2 c 2 ) w(x,p ± ), 
f(x,p) = Q(p°) 5{p 2 — m 2 c 2 ) w(x,p ± ), 



(4.38) 



where Q(p°) is the unit step function. The mean- field kinetic equations for 
these functions can be derived from Eqs. (4.30). We will give only the deriva- 
tion of the equation for f(x,p) since the equation for f(x,p) is obtained anal- 
ogously. 
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Multiplying the first of Eqs. (4.30) by 0(p°) 5(p 2 - m 2 c 2 ) and again using the 
fact that on the mass-shell = p^ j mc, we have 

9(p°) 5(p 2 - m 2 c 2 ) j? fa - e - ^V^) w = 0. (4.39) 



The transverse gradient in the momentum space, V£, can be represented as 



V^-n-A (4.40) 



where d£ = g uX d/dp x and p (| = p • n. Since u> does not depend on p p the 
operator in Eq. (4.39) may be replaced by d v p . Finally, using the relations 



d5(p 2 ~m 2 c 2 ' 



Q(p°)5(p 2 -m 2 c 2 ) = 20(p u )p 



dp 2 



and p^T^p" = 0, Eq. (4.39) takes the form 



p» (d^-J^di) f(x,p) = 0. 



(4.41) 



The analogous covariant kinetic equation for antiparticles reads 

p" + ^ 9;) /(x, p) = 0. (4.42) 

Formally, Eq. (4.41) coincides with the well-known relativistic kinetic equation 
for charged particles in a prescribed electromagnetic field (see, e.g., [10]). 

Here one comment is in order. We see that the invariant distribution functions 
f(x,p) and f(x,p) satisfy kinetic equations (4.41) and (4.42) which do not give 
any indication of the family of hyperplanes a used in the derivation of these 
equations. Note, however, that a unique solution of these equations exists only 
if f(x,p) and f(x,p) are specified on some space-like surface a in Minkowski 
space. To formulate this "initial condition" , we have to recall Eqs. (4.38) which 
relate the invariant distribution functions to the functions w and w defined on 
the family of hyperplanes o n T . Lorentz invariance of the theory manifests itself 
by associating f(x,p) and f(x,p) with an arbitrary family of hyperplanes in 
order to fix the "initial condition" . 
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5 Summary and outlook 



Based on the general density matrix approach to QED plasmas [1], we have 
derived kinetic equations for the fermionic subsystem in the mean-field approx- 
imation. The general mean-field expression given by Eq. (3.5) is a covariant 
generalization of the kinetic equation derived previously by Bialynicki-Birula 
et al. [5]. Their result is reproduced in the "instant frame", where the nor- 
mal vector is given by n M = (1, 0, 0, 0). The covariant structure of Eq. (3.5) is 
particularly convenient for the spinor decomposition which allows to separate 
kinetic processes in different channels. Another advantage of this equation is 
that it can be used to perform further approximations in a covariant form. 
For instance, the quasi-classical limit of the mean-field kinetic equation is 
presented. 

Applications of relativistic mean-field theories have been discussed in different 
contexts. For instance, in heavy-ion collisions a relativistic kinetic equation in 
an "instant frame" is solved using the relativistic Landau- Vlasov method [12]. 
Present ultra-relativistic heavy-ion collisions demand a consistent relativistic 
approach to nonequilibrium evolution [7]. Laser-plasma interactions are most 
often treated within particle in cell (PIC) simulations [13], which follow from 
classical mean-field approximations to the relativistic kinetic equation. 

Spectral information is not contained in the description presented here. In [14,15], 
for instance, the relation between the one-time and two-time Wigner function 
in the instant frame is discussed. Applying an energy moment expansion of the 
two-time Wigner function, the one-time Wigner function is given by the lowest 
moment, whereas the spectral information is contained in higher moments. 

There are different ways to go beyond the approximations presented in this 
paper. Quantum corrections to the quasi-classical kinetic equation can be 
taken into account by expanding Eqs. (3.24) - (3.31) in terms of Planck's 
constant. This implies to consider non-local fluctuations in the plasma at a 
length scale less than the de Broglie wavelength. Quantum effects associated 
with particle-antiparticle coherence, like for instance pair production caused 
by strong fields [7,8], can give significant corrections to the Wigner function 
in the treatment of QED plasmas under extreme conditions. Pair creation in 
current laser-plasma experiments is realized through a bremsstrahlung con- 
version of MeV electrons into MeV photons [11]. For the description of such 
effects the photon kinetics has to be included self-consistently into the picture 
presented here. 

Furthermore, improving the mean-field approximation, one can consider col- 
lisions in the plasma. This can be done systematically by expanding the col- 
lision terms for the photons and electrons (compare Eq. (1:5.34) and (1:5.35)) 
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in terms of the fine structure constant a. First order effects related to emision 
and absorption of photons is subject of forthcoming studies. 



Appendix A 

Kinetic equation for the matrix Wigner function in the "instant frame" 

We consider Eq. (3.5) in the "instant frame", where n M = (1,0,0,0) and 
r = x° = t (c = h = 1). Introducing the usual space-time notation x^ = (t, r), 
we find 

a£ = (0,r), x ±ll = (0,-r), tf = (0,p), P ±li = (0,-p). 

The transverse four-gradients with respect to space-time and the momentum 
variables, V„ and VJf, are written in the "instant frame" as 

v M = (o,v), y = (o,-v), v PM = (o,a p ), V£ = (o,-0 P ), 

where V = d/dr and d p = d/dp. For definiteness, Cartesian components of 
all three-dimensional vectors and gradients will be written with upper Latin 
indices running from 1 to 3. For instance, V* = d/dr l and d l p = d/dp 1 . Sum- 
mation over repeated Latin indices is implied. 

The total electric and magnetic fields, S and B, are defined in the "instant 
frame" as 

£ = -^p - VA°, B = V x A, (A.l) 
at 

or, 

= B i = £ i3k v j A k^ ( A2 ) 

where e^ k is the three-dimensional antisymmetric symbol with e 123 = I. Note 
also that, in our notation, the non-zero components of the tensor (3.9) are 
now given by Fl = F Uj = - (V i A j - V j A l ). 

Further it is easy to verify that the operators (3.6) - (3.8) can be written as 
D T = D t , D ±M = (0,D), P ±m = (0,-P), (A.3) 

where 
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1/2 

Dt = §i + C I ds£ ^ r + isd p)- d P ^ ( A - 4 ) 

-1/2 
1/2 

D = V + e J dsB(t,r + isd p )xd p , (A.5) 

-1/2 
1/2 



P 



p-ie J sdsB(r + isdp} x d p . (A.6) 

-1/2 

Finally, in the "instant frame" we have 7y = 7 = j3, so that the matrices S* M 
[see Eq. (3.4)] can be written in terms of the Dirac a-matrices as 

S li = {0,-ia). (A.7) 



Putting expressions (A. 3) and (A.7) into Eq. (3.5), it is convenient to rewrite 
this equation for the modified Wigner function [5] 

W = (A.8) 



which implies that the fermionic density operator is defined as [cf. Eq. (2.2)] 

Paa>{ry) = -\mr)^ a ,{r')\. 
Then, in terms of W, Eq. (3.5) becomes 

D t W = -im \/3, w]-]-D- {a, w\ - iP ■ [a, w] . (A.9) 

This is the mean-field kinetic equation derived by Bialynicki-Birula et al. [5]. 

Quasi- classical kinetic equations in the "instant frame" 



We now aim to show that, in the "instant frame", the kinetic equations (4.25) 
take the form (4.26). First we note that in this frame the transverse four- 
velocity (4.23) is written as v± = (0,v), where v = c 2 p/e(p) is the three- 



dimensional velocity vector of a particle with the energy e(p) = c\jp 2 + m 2 c 2 . 
Then, using the above expressions for the transverse four-gradients in the 
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"instant frame" , we find that 



d 

y _ _f if — _£ . 

"SVJ = = - (« x B) 



<9p 



(A.10) 



Finally, in the "instant frame" we have the obvious relation 

8 vt „ 1/9 \ 
tt- + — V, = - U- + « • V . 



dr 



c \dt 



(A.H) 



Insertion of expressions (A. 10) and (A. 11) into (4.25) leads to the kinetic 
equations (4.26). 



Appendix B 



Here we give some basic relations, which are used to transform the matrix 
kinetic equation (3.5) into Eqs. (3.16) - (3.20). We follow the notation of [16]. 

The totally antisymmetric Levi-Civita tensor e^ -! 3 is defined through even 
and odd permutations of with 



£ 0123 _ _£ 0123 _ \ 



(A.l) 



The relation between co- and contravariant components follows from the met- 
ric g^ = diag(l, -1, -1, -1) 

From the Dirac algebra (see e.g. [16]) we can calculate the commutator and an- 
ticommutator relations appearing in the different spinor channels of Eq. (3.5) 



<V>75 



0, 



c V'7s7 / / 



2 ^ (^'TsT/x - > 



— — 2i (g^a^, + g uu >o^i — g^ u i(y un> ~ fiV ^') > 



{ V ' 7 } = 2 V ' { V ' V } = 2^^757" > 



21 





} = 2 \9w'9w> 



{<V> 757,1'} = 

9 l xu'9v i j) I + ^ £ nuti'v' 75- 
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